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Abstract 

A variety of astronomical phenomena appear to not satisfy the ergodic 
hypothesis in the relevant stationary state, if any. As such, there is no reason 
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for expecting the applicability of Boltzmann-Gibbs (BG) statistical mechan- 
ics. Some of these phenomena appear to follow, instead, nonextensive sta- 
tistical mechanics. In the same manner that the BG formalism is based on 
the entropy Sbg = —kJ2iPi^ n Pi, the nonextensive one is based on the form 
S q = k(l—J2iPi)/(q—l) (with S\ = Sbg)- The stationary states of the former 
are characterized by an exponential dependence on the energy, whereas those 
of the latter are characterized by an (asymptotic) power-law. A brief review 
of this theory is given here, as well as of some of its applications, such as the 
solar neutrino problem, polytropic self-gravitating systems, galactic peculiar 
velocities, cosmic rays and some cosmological aspects. In addition to these, 
an analogy with the Keplerian elliptic orbits versus the Ptolemaic epicycles 
is developed, where we show that optimizing S q with a few constraints is 
equivalent to optimizing Sbg with an infinite number of constraints. 



I. INTRODUCTION 



Connections between dynamics and thermodynamics are far from being completely clari- 
fied. For instance, long-range interactions are expected to substantially modify various usual 
thermodynamical properties. E. Fermi addressed such question with the following words [1]: 

The entropy of a system composed of several parts is very often equal to the sum of the 
entropies of all the parts. This is true if the energy of the system is the sum of the energies 
of all the parts and if the work performed by the system during a transformation is equal to 
the sum of the amounts of work performed by all the parts. Notice that these conditions are 
not quite obvious and that in some cases they may not be fulfilled. Thus, for example, in the 
case of a system composed of two homogeneous substances, it will be possible to express the 
energy as the sum of the energies of the two substances only if we can neglect the surface 
energy of the two substances where they are in contact. The surface energy can generally be 
neglected only if the two substances are not very finely subdivided; otherwise, it can play a 
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considerable role. 

In those of L. Tisza [2]: 

The situation is different for the additivity postulate P a2, the validity of which cannot 
be inferred from general principles. We have to require that the interaction energy between 
thermodynamic systems be negligible. This assumption is closely related to the homogeneity 
postulate P dl. From the molecular point of view, additivity and homogeneity can be expected 
to be reasonable approximations for systems containing many particles, provided that the 
intramolecular forces have a short range character. 

Finally, in those of P.T. Landsberg [3]: 

The presence of long-range forces causes important amendments to thermodynamics, 
some of which are not fully investigated as yet. 

In recent papers, also E.G.D. Cohen [4] and M. Baranger [5] have addressed this ques- 
tion. Indeed, a significant amount of systems, e.g., turbulent fluids ( [6,7] and references 
therein), electron-positron annihilation [8], quark-gluon plasma [9], economics [10-12], mo- 
tion of Hydra viridissima [13], kinetic theory [14], classical chaos [15], quantum chaos [16], 
quantum entanglement [17], anomalous diffusion [18], long-range-interacting many-body 
classical Hamiltonian systems ( [19] and references therein), internet dynamics [20], and 
others, are known nowadays which in no trivial way accomodate within BG statistical me- 
chanical concepts. Systems like these have been handled with the functions and concepts 
which naturally emerge within nonextensive statistical mechanics [21-23]. 

The basic building block of nonextensive statistical mechanics is the nonextensive entropy 



The entropic index q characterizes the statistics we are dealing with; q = 1 recovers the 
usual BG expression 



[21] 




(q e 11). 
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We may think of q as a biasing parameter: q < 1 privileges rare events, while q > 1 privileges 
common events. Indeed, p < 1 raised to a power g < 1 yields a value larger than p, and the 
relative increase p 9 /P = p 9-1 is a decreasing function of p, i.e., values of p closer to (rare 
events) are benefited. Correspondingly, for q > 1, values of p closer to 1 (common events) 
are privileged. Therefore, the BG theory (i.e., q — 1) is the unbiased statistics. A concrete 
consequence of this is that the BG formalism yields exponential equilibrium distributions 
(and time behavior of typical relaxation functions), whereas nonextensive statistics yields 
(asymptotic) power-law distributions (and relaxation functions). Since the BG exponential 
is recovered as a limiting case, we are talking of a generalization, not an alternative. 

To obtain the probability distribution associated with the relevant stationary state (ther- 
mal equilibrium or metaequilibrium) of our system we must optimize the entropic form (1) 
under the following constraints [21,22]: 

£ft = 1 > (3) 

i 

and 

^iPi Ei = tj u\ 

where {Ei} is the set of eigenvalues of the Hamiltonian (with specific boundary conditions), 
and U q is a fixed and finite number. This optimization yields 

„ _ [I - (I - qWEj - U q )]^ 

Pi — ^ 5 K°) 

where 

Z q = £[1 - (1 - q)Pq(E 3 - U q )}^ , (6) 



and 



Erf 



h = ( 7 ) 



P being the optimization Lagrange parameter associated with the generalized internal energy 
U q . Equation (5) can be rewritten as 



Pi oc [1 - (1 - qtfEi] 1 " 1 -* = e~/ E * , (8) 

where f3' is a renormalized inverse "temperature" , and the q- exponential function is defined 
as e x = [1+(1 — q)x] 1 ^ 1 ~ q ^ = 1/[1 — (?— l)^] 1 ^ 9-1 ^ (with ef = e x ). This function replaces, in a 
vast number of relations and phenomena, the usual BG factor. In particular, the ubiquitous 
Gaussian distribution oc e~ ax2 becomes generalized into the distribution oc e~ aqx2 = 1/[1 + 
(q - l)a g z 2 ] 1/(9-1) (fat-tailed if q > 1, and with a cutoff if q < 1). 

II. ANALOGY WITH THE PTOLEMY-KEPLER PROBLEM 

Ancient Greeks believed, through philosophical arguments, that the orbit of any celestial 
body ought to be a "perfect" geometrical figure, namely a circle. This belief was very 
consistent with the observations of the motion of the stars. Not so, however, with the 
observations of the motion of the planets (a word whose etymology precisely is wanderer). 
This nontrivial problem was attacked by Ptolemy. He introduced, with sensible success at 
that time, the idea that the motion of a planet would be that of a circle rotating around 
another circle, in turn rotating around another circle, and so on. The idea of such epicycles 
was recursively used by most astronomers along the centuries over large regions of the 
civilized world [24]. Before Kepler's time, some astronomers had used up to several dozens(!) 
of epicycles one on top of the other, in order to increase the precision of their calculations. 
Then Kepler arrived and proposed that the orbit of a planet around the sun would be that of 
a geometrical figure which generalizes the circle, namely that of an ellipse. This figure needs, 
for being characterized, one new parameter, namely the eccentricity e, in addition to say its 
largest diameter. The circle is thus recovered as the particular case e = 0. In very few years, 
practically all the astronomers of Europe abandoned the cumbersome calculations with the 
Ptolemaic epicycles, and adopted the Keplerian ellipses. Ptolemy idea was nevertheless not 
wrong at all! Indeed, the planetary elliptic motion can be described in terms of infinite 
series, whose terms can be identified with the Ptolemaic epicycles. In other words, the 
elliptic keplerian motion can be regarded as resulting from an "infinite series of epicycles" . 
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A very lucid discussion of the correspondence between (i) the series expansion of the elliptic 
motion and (ii) the Ptolemaic geometric construction, was provided by Hoyle [25]. The 
alluded series for the planetary elliptic motion can be found in any text-book of Celestial 
Mechanics [26] , and are given (in polar coordinates, as a function of time t) by 

5 

9 = nt + 2esinnt + - e 2 sin2nt + ... (9) 
4 

and 

r 1 

- = l - ecosnt + - e 2 (1 - cos2nt) + ... (10) 
d Zi 

where a is the semi-major diameter, and n = 2n/P, P being the so called sidereal period. 
If the orbit is a circle, i.e., e = 0, then each polar coordinate has only one term in its 
expansion. For any value of e ^ 0, an infinite number of terms must be used if one wants 
the exact Keplerian answer. This sudden change in the number of terms of course reflects 
the equally sudden change of symmetry when one goes from the circle to the ellipse. In 
other words, if one wants to stick to circles, that is perfectly possible and correct, although 
the price to be paid is that one has to consider an infinite number of circles. The other, 
simpler, possibility, is of course to use only one ellipse. One has however to know, through 
any accessible procedure, the value of e for that particular planet. 

Let us go back now to the problem of the use and status of the nonextensive entropy 
S q as given in Eq. (1). How universal is the BG entropy? How S q fundamentally relates 
to it? What is the status of the entropic index q, which characterizes universality classes of 
nonextensivity? (The most famous one being of course the q — 1 universality class). Let us 
address this epistemologically interesting point through the canonical ensemble associated 
with conservative Hamiltonian systems. 

Let us first develop the q-logarithm function 



In x = x ' x = ^_ 1 = \ nx _ i n 2 x + ln s x+ (n) 

1 — q 1 — q 2 6 



1-9 _1 e (l-q)lnx_ 1 5-1,2 (<? ~ X ) 
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as well as 



p 9-l = e fe-l)ln W = 1 + ( g _ 1 ) lnft+ (9 y! ln 2 p . + ( 12 ) 

The functional to be optimized for the canonical ensemble is 

- + «& + ^ = + « + > (13) 

where a and /3 are Lagrange parameters, and (...) = J2iPi(---)- We use now the series (11) 
and (12) inside Eq. (13), and obtain 

*■({*}) = Sj ^^ - ^Vft> - ^^fln'ft) + ... + a £p< 

+/3<£ i ) + /% - 1)[(^ lnp,) - (^)(lnp,)] 
+/5^-y^[(^ln 2 ^) - (^)(ln 2 Pl ) + 2(^)(ln Pl ) 2 - 2(E t \n Pi ) (\np t )] + ... (14) 

= *i({ft» 

" " V^) - ii^ll!(ln 3 ^) + ... + l)[(^ln ft > - (^)(lnp,)] 



o 1)2 [(E i ln 2 p i ) - (^)(ln 2 Pl ) + 2(^)(lnp,) 2 - 2(E l \n Pi ) (\n Pi )] + ... (15) 



2 *■ 1 Q 

2 

This optimization is the usual one within BG statistical mechanics if q — 1, having a /eiu 
constraints. It suddenly becomes a quite complex one, with an infinite number of constraints, 
if q 7^ 1. Let us further analyze the generic case. 

Let us consider the following generic optimization problem on the basis of the Sbg 
entropy. The constraints to be used are Eq. (3), as well as the following ones: 

(ln k p) = C k (A; = 2,3,4,...), (16) 

where {Ck} are fixed finite quantities, to which we shall associate the Lagrange parameters 
{7a,.}. In addition to these constraints, we also have the following ones: 

(Ei) = D 1 , (17) 
(Eihipi) - (Ei)(]n Pi ) = ((Ei - (Ei))Qn Pi - (lnp,))) = D 2 , (18) 



(EM 2 Pr) - (Ei)(ln 2 Pi ) + 2(E l )(\np l ) 2 - 2(E i lnp i )(]np i ) = D 3 , (19) 

7 



etc, where Di, D 2 , D 3 , Di, ... are fixed finite quantities, to which we shall associate the 
Lagrange parameters 5i,5 2 ,5 3 , ...,Si, .... The optimization of Sbg with all these constraints 
provides an optimizing distribution noted p(Ef, a, {7fc}, {Si}). We shall next consider a 
particular case of this distribution, the one where we impose 

(q - l)*" 1 

lk = - {q J (A; = 2,3,4,...), (20) 

61 = P, (21) 
6 2 = (3(q-l), (22) 

* = />^, (23) 

etc. Let us stress that taking a particular case of p(Ei, a, {7^}, {5;}) is perfectly legitimate. 
[Indeed, as an illustration, let us remind that in the usual BG statistical mechanics, the 
distribution corresponding to the microcanonical ensemble can be viewed as the (3 = 
particular case of the distribution associated with the canonical ensemble]. Consequently 
the stationary distribution becomes a function only of (Ei,(3,q). We may alternatively 
consider that we have defined a specific new ensemble, whose distribution is fully known 
once we know (3 and q. The inverse "temperature" (3 depends on the physical thermostat 
with which the system is in contact, and q depends on the microscopic dynamics of the 
specific physical system. The symmetry change that has occurred from q = 1 to q 7^ 1 
hopefully corresponds to the fact that ergodicity (homogeneity in phase space) has been 
lost, and occupancy is now related to a presumably (multi)fractal or hierarchical geometry 
in phase space. The fact that only one new parameter emerges (namely q), and not the 
infinite number that would correspond to the entire multifractal function f(a), should be 
presumably related to the expectation that the thermodynamical properties would basically 
depend only on one number, namely the Hausdorff dimension of the system (see [27] for a 
connection between q and the Hausdorff dimension for the simple case of the logistic-like 



maps). Although further studies are necessary to transparently establish the connection of 
the geometry in phase space with q, we have shown herein that indeed the optimization of 
S q with a few constraints is equivalent to the optimization of Sbg with an infinite number of 
them. Therefore, it is in principle possible and correct to use Sbg for such complex systems. 
The price to be paid is the fact that an infinite number of constraints must be, in one way 
or another, taken into account. The alternative is to use S q with just a few constraints, 
but one needs to know the value for q to be used. The analogy with the Ptolemy-Kepler 
problem, where e must be known, is therefore, in our opinion, quite striking. 

III. APPLICATIONS 

Let us now briefly review five recent applications of the ideas associated with nonextensive 
statistical mechanics to phenomena in astronomy and astrophysics, namely the solar neutrino 
deficit [28], self-gravitating polytropic systems [29-31], peculiar velocities of galaxy clusters 
[37] , the flux of cosmic rays [38] , and some cosmological effects [39] . 
Solar neutrino problem: 

The solar plasma is believed to produce large amounts of neutrinos through a variety of 
mechanisms (e.g., the proton-proton chain). The calculation done using the so called Solar 
Standard Model (SSM) results in a neutrino flux over the Earth, which is roughly the double 
of what is measured. This is sometimes referred to as the neutrino problem or the neutrino 
enigma. There is by no means proof that this neutrino flux defect is due to a single cause. It 
has recently been verified that neutrino oscillations do seem to exist, which would account for 
part of the deficit. But it is not at all clear that it would account for the entire discrepancy. 
Quarati and collaborators [28] argue that part of it - even, perhaps, an appreciable part of 
it - could be due to the fact that BG statistical mechanics is used within the SSM. The solar 
plasma involves turbulence, long-range interactions, possibly long-range memory processes, 
all of them phenomena that could easily defy the applicability of the BG formalism. Then 
they show [28] in great detail how the modification of the "tail" of the energy distribution 
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could considerably modify the neutrino flux to be expected. Consequently, small departures 
from q — 1 (e.g., \q — 1| of the order of 0.1 or even less) would be enough to produce as much 
as 50% difference in the flux. This is due to the fact that most of the neutrino flux occurs at 
what is called the Gamow peak. This peak occurs at energies quite above the temperature 
(say 10 times larger), i.e., at energies in the tail of the distribution. 
Polytropic Equilibrium Solutions to the Vlasov-Poisson Equations: 

The first physical application of the non-extensive thermostatistical formalism was re- 
lated to the study of maximum entropy solutions to the Vlasov-Poisson equations describing 
self gravitating iV-body systems like galaxies [29,30]. The maximization of the standard 
Boltzmann-Gibbs entropy under the constraints imposed by mass and energy conservation 
lead to the isothermal sphere distribution, which has infinite mass and energy [32]. In [29,30], 
it was shown that the extremalization of the non extensive g-entropy under the same con- 
straints leads to the stellar polytropic sphere distributions which, for a certain range of the 
q parameter, are endowed with finite mass and energy, as physically expected. This consti- 
tuted the first clue suggesting that the generalized thermostatistical formalism based on S q 
may be of some relevance for the study of systems exhibiting non extensive thermodynam- 
ical properties due to long range interactions. The possible usefulness of non-logarithmic 
entropic measures in the study of stellar systems had also been suggested in [33]. 

Stellar polytropic sphere distributions are of the form 

/(x, v) = /(e) = A($ - e)"- 3 / 2 e < $ (24) 

e > $ , 

where 

e = iv 2 + $(x), (25) 

is the total energy (per unit mass) of an individual star, and A, $ , and n (usually called 
polytropic index) are constants. The quantity /(x, w)d 3 xd 3 v denotes the number of stars with 
position and velocity vectors respectively within the elements d 3 x and d 3 v in position and 
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velocity spaces. The polytropic distribution (24) exhibits, after an appropriate identification 
of the relevant parameters, the g-MaxEnt form. The entropic parameter q can be expressed 
in terms of the index n by identifying n — 3/2 with 1/(1 — q) (see [30] for details; notice that 
here we are identifying / with p, and not with the escort distribution oc p q ) . We obtain 

2n - 5 

" = 2^3- (26) 

The limit n — > oo (hence q — 1) recovers the isothermal sphere case; for n < 5 (hence 
q < 5/7), finiteness for mass and energy naturally emerges. Polytropic distributions consti- 
tute the simplest, physically plausible models for self-gravitating stellar systems [32]. Alas, 
these models do not provide an accurate description of the observational data associated with 
real galaxies. In spite of this, the connection between the S q entropy and stellar polytropes 
is of considerable interest. Besides the (in our opinion) notable fact that, for a special range 
of values of q, non-extensive thermostatistics leads to finite stellar systems, the established 
connection between the S q entropy and stellar polytropic distributions is interesting for other 
reasons. Polytropic distributions arise in a very natural way from the theoretical study of self 
gravitating systems. The investigation of their properties has been of constant interest for 
theoretical astrophysicists during the last one hundred years [32]. Polytropic distributions 
are still the focus of an intense research activity [31], and the study of their basic proper- 
ties constitutes a part of the standard syllabus of astrophysics students. Now, polytropic 
distributions happen to exhibit the form of q-MaxEnt distributions, that is, they constitute 
distribution functions in (x, v) space that maximize the non-extensive functional S q under 
the natural constraints imposed by the conservation of mass and energy [29] (other con- 
straints associated with other conserved quantities may be incorporated). Let us recall that 
the original path leading to the S q entropic form was not motivated by self-gravitating sys- 
tems, nor was it motivated by any other specific system. The entropic form S q was proposed 
by recourse to very general arguments dealing with the consideration of (i) entropic forms 
incorporating power law structures (inspired on multifractals) and reducing to the standard 
logarithmic measure in an appropriate limit and (ii) the basic properties a functional of 
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the probabilities should have in order to represent a physically sensible entropic measure 
[21]. The g-entropy is a very natural and, to a certain extent, unique generalization of the 
Boltzmann-Gibbs-Shannon measure. Taking this into account, it is remarkable that the 
extremalization of the g-measure leads to a family of distribution functions of considerable 
importance in theoretical astrophysics. In a sense, we might say that astrophysicists, when 
studying newtonian self-gravitating systems, have been working with g-MaxEnt distribution 
functions for a hundred years without being aware of it. The discovery of the connection 
between S q and stellar polytropic distributions stimulated the application of nonextensive 
statistical mechanics to the study of other systems with long range interactions [34,35]. In 
particular, the analysis by Boghosian [34] of metastable states in pure electron plasmas lead 
to the identification of the first maximum g-entropy distribution observed in the laboratory 
[34,35]. In point of fact, the mathematical formalism used by Boghosian to describe the 
electron plasma is closely related to one associated to the stellar polytropic distributions 
[34]. 

Peculiar velocities of galaxy clusters 

The COBE (Cosmic Background Explorer) satellite measured the peculiar velocities 
(difference of velocity with regard to the average expansion of the universe) of some clusters 
of spiral galaxies. A distribution was found which exhibits a cutoff around 500 Km/s. The 
Princeton astrophysical group [36] analyzed the distribution of velocities within four different 
cosmological models (cold matter, hot matter, premieval barionic). None of those attempts 
succeeded in reproducing the observed cutoff, although each of those models involved several 
free parameters (that were fixed through a variety of arguments). By assuming within 
nonextensive statistical mechanics, an extremely - almost outrageously - simplified model, 
namely an ideal classical gas!, the empirical velocity distribution was quite satisfactorily 
matched [37]. Only two fitting parameters were used, namely the scale of velocities and 
q ~ 0.23. In spite of the extreme simplicity of the model, the fact that the statistics 
was allowed to change proved its high efficiency The possible primacy of statistics over 
models is not new in statistical physics in fact. Even if we consider the ideal gas hypothesis 
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for a quantum gas, we obtain acceptably good first approximations for phenomena such 
as superconductivity or suprafluidity (using Bose-Einstein statistics), or such as metallic 
conductivity of electrons (using Fermi-Dirac statistics). Of course, more realistic models 
must include interactions, but an acceptable first approach can indeed be done already on 
the basis of an appropriate statistics (even if the model is oversimplified). 
Cosmic rays: 

The flux of cosmic rays arriving on Earth is a quantity whose measured range is among 
the widest experimentally known (33 decades in fact). This distribution refers to a range of 
energies which also is impressive (13 decades). This distribution is very far from exponential, 
as can be verified on Fig 1. This basically indicates that no BG thermal equilibrium is 
achieved, but some other (either stationary, or relatively slow varying) state, characterized 
in fact by a power law. If the distribution is analyzed with more detail, one verifies that 
two, and not one, power-law regimes are involved, separated by what is called the "knee" 
(slightly below 10 16 ev). At very high energies, the power-law seems to be interrupted by 
what is called the "ankle" (close to 10 19 ev). One may guess that, along such wide ranges 
(of both fluxes and energies), a variety of complex intra- and inter-galactic phenomena 
are involved, related to both the sources of the cosmic rays as well the media they cross 
before arriving on Earth. However, from a phenomenological viewpoint, the overall results 
amount to something quite simple. Indeed, by solving a simple differential equation, a quite 
remarkable agreement is obtained [38]. This differential equation is the following one: 

|| = -b'pt - ftp? . (27) 

This differential equation has remarkable particular cases. The most famous one is (</, q) = 
(1,2), since it precisely corresponds to the differential equation which enabled Planck, in 
his October 1900 paper, to (essentially) guess the black-body radiation distribution, thus 
opening (together with his December 1900 paper) the road to quantum mechanics. The more 
general case q' = 1 and arbitrary q is a simple particular instance of the Bernoulli equation, 
and, as such, has a simple explicit solution. This solution has proved its efficiency in a variety 
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of problems, including in generalizing the Zipf- Mandelbrot law for quantitative linguistics 
(see Montemurro's article in [23] for a review). Finally, the generic case q > q' > 1 also 
has an explicit solution (though not particularly simple, but in terms of two hypergeometric 
functions) and produces, taking also into account the ultra-relativistic ideal gas density of 
states, the above mentioned quite good agreement with the observed fluxes. Indeed, if we 
assume < b' « b and q' < q, the distribution makes a crossover from a power-law 
characterized by q at low energies to a power-law characterized by q' a t high energies, which 
is exactly what the cosmic rays exhibit to a quite good approximation: See Fig. 1. 
Cosmology: 

Nonextensive statistical mechanics has also been applied to a variety of cosmological 
and general relativity problems including the cosmic background radiation in a Robertson- 
Walker universe, the dynamics of inflationary cosmologies, the universal density profile of 
dark halos, early universe phenomena (e.g., the primordial 4 He formation), among others 
[39]. 

IV. CONCLUSIONS 

In outer space physics and astrophysics, there is a considerable amount of anomalous 
phenomena that require a thermostatistical treatment which exceeds the usual capabilities 
of Boltzmann-Gibbs statistical mechanics. This fact is due to the relevance of gravitational 
forces (which are long-ranged), as well as to a variety of dynamical nonlinear dynamical 
aspects. Some of these phenomena appear to be tractable within nonextensive statistical 
mechanics, and we have illustrated this with a few typical examples. 

In addition to this, we have argued that, in analogy with the Ptolemy-Kepler problem 
concerning the orbits of the planets, the optimization of the nonextensive entropy S q with 
a few constraints can be equivalently seen as the optimization of the BG entropy with an 
infinite number of (ultimately related) constraints. The Keplerian orbit is equivalent to 
infinite Ptolemaic epicycles, and its series expansion in powers of the eccentricity e abruptly 
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reduces to a single cercle when e precisely vanishes. Quite analogously, the variational 
problem associated with S q for the canonical distribution is equivalent to considering Sbg 
with infinite (related) constraints, and its series expansion in powers of q—1 abruptly reduces 
to the Sbg variational problem with a single constraint (besides the trivial one associated 
with the norm) when q — 1 precisely vanishes 

Last but not least, the statements by Fermi and by Tisza that we quoted in the Intro- 
duction constitute a reply which precisely points one of the various inadvertences contained 
in Nauenberg's recent criticism of nonextensive statistical mechanics [40]. The nonextensive 
thermostatistical formalism addresses, among others, conservative Hamiltonian many-body 
systems including long-range interactions. For such systems, the "surface energy" (using 
Fermi's expression), or the "interaction energy" (using Tisza's expression), is as relevant as 
the energies of the parts, and does not become negligible in the thermodynamic limit, in 
contrast with what happens for short-range interacting systems. Nauenberg violates pre- 
cisely this by imposing his Eq. (7). Moreover, a few lines earlier (his Eq. (5)), he imposes 
the factorization of the probabilities. This property is inconsistent with his own Eq. (7). 
Indeed, the metaequibrium distribution of energies is given by the g-exponential function 
e x q = [1 + (1 - g)x] 1/(1_9) , and, unless q = 1, generically e x q +y ^ e x q e\ = e *+y+^-i)*y . This 
type of inconsistency (for finite systems) has already been pointed out by several authors, 
starting with [41]. It is not impossible that this specific inconsistency disappears in the limit 
iV — > oo, if taken before the t — > oo limit, but this is highly nontrivial and remains to be 
proved. 
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FIG. 1. Flux of cosmic rays as a function of their energies. See [38] for the origin of the exper- 
imental data and theoretical details. The blue dashed curve corresponds to the BG distribution, 
and the red continuous curve to a crossover between a q = 1.225 distribution (before the knee) and 
a q' = 1.185 one (after the knee). 
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